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ABSTRACT 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 



The problem of detecting an enemy submarine whose possible position was 
revealed by the hit of a torpedo is known as the “Flaming Datum” problem. All previous 
studies devoted to this theme make unrealistic assumptions about the speed of the 
escaping target when dealing with a diesel-electric submarine. In this kind of submarine 
the constraint imposed by the remaining charge of its batteries determines that its 
behavior is essentially conservative in how fast it should escape. 

The objective of this thesis is to explore the idea of varying speed in the flaming 
datum problem. Two different approaches are considered. An analytical model is 
developed based on the relationship among some of the physical factors that could 
determine or constrain the behavior of a diesel submarine while escaping from the area of 
the flaming datum. The second approach considers a discrete event simulation using the 
Java-based Simkit package. Data analysis is used to determine a possible fit for the 
simulation results. Several tactics are explored to determine their effects on detection 
probability. 
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EXECUTUVE SUMMARY 



The problem of detecting an enemy submarine whose possible position has been 
revealed by the hit of a torpedo on a surface vessel is known as the “Flaming Datum” 
problem. Since World War II, several studies have been devoted to the issue where a 
submarine tries to escape from the datum area while being hunted by a helicopter 
equipped with variable depth sonar (VDS) and/or antisubmarine weapons. All these 
studies make assumptions about the speed of the escaping target that are unrealistic when 
dealing with a diesel-electric submarine. The constraint imposed by the remaining 
charge of its battery determines that the behavior of this kind of submarine is essentially 
conservative in how fast it should escape. Even though the submarine is able to maintain 
a constant speed for long periods of time, this is only valid for very low speeds. Low 
speeds are not tactically acceptable for a submarine trying to evade a searcher, especially 
in the first minutes after the attack when operating in the vicinity of the Datum. 

According to opinions and experiences collected by the author, a diesel-electric 
submarine will tend to escape, as a general rule, as soon as possible from the zone where 
an attack recently occurred. This tactic is motivated by the efforts of the enemy to find the 
submarine, but it is only sustainable for short periods of time depending on the 
characteristics of the battery, and its remaining charge. 

The objective of this thesis is to develop and test a model for the Flaming Datum 
problem with a submarine trying to evade at varying speeds. First, the author develops 
an analytical model based on the existing relationship among some of the physical factors 
that determine or constrain the speed and behavior of a diesel submarine trying to escape 
from the area of the flaming datum. This model predicts the probability of detecting the 
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submarine as a function of the helicopter’s delays in arriving at the datum area. To 
obtain a more accurate result for the same problem a discrete event simulation using the 
Java based Simkit package is developed and validated as believable. This simulation is 
also used to test some variations where either the submarine or the helicopter tries a 
specific tactic. 

Even though, both the analytical model and the simulation show some 
weaknesses that may affect the credibility of the results, the possible sources of 
inaccuracy are identified and explained. It turns out that the submarine speed should be a 
strong function of time, rather than constant as usually supposed. 

This study represents a more realistic method of posing the Flaming Datum 
problem when dealing with diesel-electric submarines that must respect an energy 
constraint. This constraint was never considered in any other study. For this reason the 
developed models may be considered a valid starting point for future works. 
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I. INTRODUCTION 



A. DESCRIPTION OF THE PROBLEM 

The problem of detecting an enemy submarine whose possible position has been 
revealed by the hit of a torpedo on a surface vessel is known as the “Flaming Datum” 
problem. Since World War II, several studies have been devoted to the issue where a 
submarine tries to escape from the datum area while being hunted by a helicopter 
equipped with variable depth sonar (VDS) and/or antisubmarine weapons. All these 
studies make assumptions about the speed of the escaping target that are unrealistic when 
dealing with a diesel-electric submarine. The constraint imposed by the remaining charge 
of its battery determines that the behavior of this kind of submarine is essentially 
conservative in how fast it should escape. Even though the submarine is able to maintain a 
constant speed for long periods of time, this is only valid for very low speeds. Low speeds 
are not tactically acceptable for a submarine trying to evade a searcher, especially in the 
first minutes after the attack, when operating in the vicinity of the Datum. 

According to opinions and experiences collected by the author, a diesel-electric 
submarine will tend to escape, as a general rule, as soon as possible from the zone where 
an attack recently occurred. This tactic is motivated by the efforts of the enemy to find the 
submarine, but it is only sustainable for short periods of time depending on the 
characteristics of the battery, and its remaining charge. 

The objective of this thesis is to develop and test a model for the Flaming Datum 
Problem with a submarine trying to evade at varying speeds. To do so, the author 
develops an analytical model based on the existing relationship among some of the 
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physical factors that determine or constrain the speed and behavior of a diesel submarine 
trying to escape from the area of the flaming datum. The model is tested with a discrete 
event simulation. 



B. BACKGROUND 



The previous studies can be divided into those that have been used to solve the 
problem as either one-dimensional or two-dimensional games. In the first group, 
Meinardi [Ref. 1] considers a multi-staged search game where the submarine hides in a 
row of Ho boxes. At each stage, the helicopter chooses a box and inspects it. If the target 
is not found, it can move on to an adjacent box or remain in the same one. Both target 
and searcher keep track of the searched boxes. The time spent to transit among the boxes 
is not taken into consideration, and thus, the applications of this game are limited for 
further analysis. 

Baston and Bostock [Ref. 2] develop a game that could be considered a 
continuous version of the flaming datum problem in one-dimensional space. The 
submarine moves in a long narrow channel and the helicopter searches for it by carrying 
a limited number of bombs. Each bomb has a lethal radius R and a limited number of 
bombs are available. These two characteristics can be considered analogous to the sonar 
detection range and the number of dips available. In this game the submarine and the 
helicopter are restricted in movement by their maximum speeds ( U and V) and 
instantaneous changes in velocity are allowed. 

The second group deals with two-dimensional cases. Danskin [Ref. 3] supposes 
that the datum is known precisely and the submarine strategy is to choose a course 0 and 
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a speed v with 0<= v<= v 0 , where Vo is the maximum speed. It holds to both throughout 
the search and this establishes the submarine's position within its speed circle [Figure 1]. 
Danskin also supposes that the submarine does not know what the helicopter is doing, so 
that it does not need to alter its speed or course. The helicopter dips are assumed to be 
inside the submarine’s speed circle and non-overlapping. Danskin’ s assumptions seem to 
be very drastic. 

Cheong [Ref. 4] develops a game where the submarine and helicopter are 
assumed to be moving at the constant speeds of U and V respectively. However, the 
submarine is allowed to execute evasive maneuvers by moving directly away from the 
last helicopter dip as well as moving away from the datum. Meanwhile the helicopter 
conducts a random search. 




Figure 1. Danskin approach. Submarine speed circle, showing a choice 
of strategy. (From Ref. 3, p. 511) 
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Washburn [Ref. 5, p 2-7] states that expected number of detections for a random 
search in an expanding disk can be modeled as: 



where V is the searcher’s speed, R is the range of its cookie cutter sensor, x is the delay 
of the helicopter arriving at the datum, znAy(t) is the expanding disk radius at time t, 
which is dependent on the submarine’s maximum speed. Equation 1 will be used in 
developing an analytical model in Chapter II. 

Washburn and Thomas [Ref. 6] model the flaming datum problem as a dynamic 
search game. They relax the rules of motion for the submarine, so that its future position 
could be anywhere inside its speed circle and not necessarily connected to the previous 
position. This assumption, in a sense, permits an unbounded target speed. In fact, the 
maximum speed is only considered to generate the radius of the area to be searched for, 
or the maximum speed circle. Even though this assumption appears to favor the 
submarine, Washbum and Thomas are of the opinion that it is of minor consequence. The 
authors compare the similarity of their results with Cheong’s, and state that this similarity 
might imply that both models are actually wrong. According to Washbum “no Flaming 
Datum Problem has ever been solved or even usefully bounded, so the required standard 
of comparison is missing” [Ref 6. P 421]. 



00 






(Equation 1) 
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c. 



APPROACH 



Our first approach is to use an analytical model to develop a plausible speed 
profile for the submarine. This model is based on the existing relationship among some of 
the physical factors that determine or constrain the speed and behavior of a diesel 
submarine trying to escape from the area of the flaming datum. Parameters and variables 
such as power, energy, speed, speed sustainability and time are employed to determine an 
expression to compute the probability of detection of the submarine as a function of time. 
The submarine’s foe is assumed to be a generic helicopter equipped with variable depth 
sonar (VDS) with a “cookie-cutter” sensor. This helicopter arrives at the flaming datum 
after some delay r and has a limited endurance. If the submarine is not detected after 
this period, it will be considered to have escaped. 

In a second and more accurate approach a Monte Carlo discrete event simulation 
using the Java-based Simkit Software package is developed. A submarine with a set of 
continuous decreasing velocity functions that are consistent with the capabilities of a 
generic submarine used in the analytical model is considered. 

A detailed description of the assumptions and strategies for the submarine and 
helicopter are given in the next chapters. 
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II. ANALYTICAL MODEL 



A. ASSUMPTIONS 



As mentioned in Chapter I, all previous studies make assumptions about the 
speed of the target that are unrealistic when dealing with a diesel-electric submarine. The 
constraint imposed on this kind of submarine by the remaining charge of the batteries 
forces its behavior to be essentially conservative when considering escape speed. We 
make the following assumptions: 

1. Submarine Assumptions 

a. The submarine will normally try to escape as quickly as possible in the 
first few minutes in order to make the FOC (farthest on circle) larger, and 
to increase the probability of escape. This behavior is justified since the 
first reaction of an attacked force will be to verify that the submarine is no 
longer in the vicinity of the datum. 

b. The battery is initially close to having a maximum charge. This is a 
reasonable assumption for a deliberate attack. 

c. The submarine is able to estimate the expected helicopter’s delay in 
arriving in the datum area. 

d. The depth factor will be ignored. 
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2. Helicopter Assumptions 

a. The position of the datum is known with no error. 

b. The searcher knows the submarine's capacities related to its speed. This 
assumption is considered valid when taking into account that open sources 
and commercial manuals are readily available. As a general rule, the 
searcher should know who the possible enemies are. 

c. The helicopter moves at a constant speed when not dipping. The times to 
dip and raise the sonar, and search in a position, are also constant. 

d. The helicopter does not assume any particular behavior for the submarine; 
therefore its dips are uniformly distributed in the submarine’s maximum 
possible distance circle (MPD). 
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B. 



THEORETICAL FOUNDATION 



Rydill and Burcher [Ref 7, p 291], state that the effective power {Pe) needed to 
“tow” the submarine at speed u{t ) in a submerged condition is given by: 



Pe - — pii{t) 



{Equation 2) 



where S REF is the wetted area of the bare hull and casing plus the wetted area of all 
appendages minus the area covered by the appendages, C TS is the total resistance 



coefficient or drag coefficient and p is the sea water density. 



Now let 
J ~ REFITS 

Thus, the energy of the submarine’s batteries required to move a submarine with a 
speed u(t) during an interval [0, t] can be expressed as: 



E{t) = J ^u{xf dx {Equation 3) 



Then, the constraint imposed by the energy of the batteries if the submarine 
wants to move for a given period [0, 7] with a speed u(t) is given by : 

E{T) < Etot {Equation 4) 



Etot is the total battery energy when the submarine begins to escape at time 0. 
This energy, according to assumption b), should be very close to maximum charge. Etot 
would be normally measured in [KWHR], but for convenience we choose units for energy 
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such that J= 1 in Equation 3. Therefore, we can compute Etot as U 3 multiplied by the 
endurance at this speed, where U is any arbitrarily chosen speed. A realistic value for 
Etot for a modem diesel-electric submarine is: 



(16.4 mphf (1.2 hr) = 5293 miles 3 hrs 2 (. Equation 5) 



C. MATHEMATICAL FORMULATION 



In Equation 1, the detection rate for a random search in an expanding disk is 
modeled as ( VW)/{7ty{t ) 2 ) . Therefore the average number of detections in the time 
interval [t, T] when the search is continuous is: 




where a = 2RV/n. 

However our helicopter searches by dipping in an expanding circle, so the value of alpha 
will depend on the position of the dips and will be changing as a function of time. 
Therefore a better expression for a(t) is: 

rate of area covered JtR 2 R 2 



a(t) = 



n 



n 







r« + £ xo'i 


V ) 




l v J 



{Equation 6) 



where C is a factor near 1 that allows computing the average distance to be flown from 
dip to dip and £is the time to make a dip. If we make C=l, the average number of 
detections in the time interval [r, 7] becomes: 
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f; 



<5 + 



R 2 

y(t ) 



-dt - 



R 1 



y( O 2 



f- 






* ,.,2 . J( 0 3 



JF 



■4- 



dt 



y(0 + 



;K'r + 



2 , >’( 0 3 



JF 



where AT = — 

8 

A logical objective function for the submarine is to minimize the expected 
number of detections during the interval [r, 7], 



PROGRAM 1 

Minimize 

1 .v(<) 

Subject to 

w(/) 3 dt < Etot and y(0) = 0 

Before time x the submarine should escape with constant speed y(r)/ r . Therefore the energy 
constraint for Program 1 is: 

r” 2 y(r ) 3 + u(t) 3 dt < Etot 
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This program can be expressed in terms of a Lagrangian unconstrained function as: 



F{y(t),u(t),t) = -^ + M0 3 
y(*Y 



(Equation 7) 



which is equal to 






K 



y( O’ + y( ' r 



+ Xu(tf 



(Equation 8) 
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Euler’s equation from the calculus of variations [Ref. 8, p 399], states that a necessary 
condition for optimality is: 



cu(t) 



(Equation 9) 



Euler also requires that 

dF 

8u(T) 



C y(T),u(T),T) = 0^u(T) = 0 



(Equation 10) 



Then 



L = 



K 



5V 



+ Xu(tf -2Xu(tf = 



K 



SV 



■ - 2Xu(t ) 3 ; r < t <T (Equation 1 1) 



and solving for u(t) 



u(t) = 





f > 




1 


K L 




2X 


y(0 2 + 

v oV ) 





T<t <T 



If we make A - — 
2/1 



and 




(Equation 12) 
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we have 



A 




( Equation 13) 



If — = Y 2 and B Ui 

B 



= v, u(t) becomes : 




(Equation 14) 



Y is an adopted constant with units of distance and v has units of speed. 

Given y(r ) , Equation 14 is a first order ordinary differential equation that can be solved 
for y(t) ; r <t <T. 

Equation 4 states that E(T ) < Etot. For this reason the submarine’s speed might 
be any value less than or equal to u(t). Therefore, the submarine’s possible locations 
could be considered uniformly distributed over the circle of radius y(t) (FOC). 
According to Washburn [Ref. 5, p. 2-10], the probability of detection as a function of 
time for a random search, given y(t), can be expressed as: 



where the random variable X is the time to detection. 

Even though Equation 14 seems not to have an analytical solution, it is possible to 
solve it numerically using the solver feature of the Excel spreadsheet. This is done by 



a(x)dx 



PD(t)=l-e 



y(x) 2 



= P(X<t);r<t<T 



(Equation 15) 
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to the constraint imposed by Equation 4, and changing the values of y(x), v and Y. 

It is unknown whether this problem possesses the convexity properties that would 
guarantee global optimality. However, when the problem was solved repeatedly using a 
wide range of different starting values for the variables y(x), v and Y, it always converged 
to the same solutions. We suspect, but can not guarantee, this is a global optimum. 

Equation 14 might be seen as a possible speed function for the submarine and 
Equation 15 is the probability of detecting it when the search is executed inside the FOC. 
However, the helicopter doesn't know what is the enemy is really doing, so it has to 
search covering every possible target’s behavior. The helicopter must assume that a 
possible objective function for the submarine is to maximize y(z) at every time t by 
using all of Etot. This assumption is expressed in Program 2, which has t as parameter for 
0 <t<T: 

PROGRAM 2 
Maximize 



Since the integrand is a convex function of u(x), the speed should be constant at u{t) for 
0 < x < t , and 



y(0 



Subject to 




tu\t) 3 = Etot-, 0 < t <T 



{Equation 16) 



therefore 



14 



u (t) = (Etot / /)' 3 ;0</<J 



{Equation 17) 



Multiplying by t 

u’ {t)t = y'(t) = E T oT U 3 t 2/3 ;0 <t<T {Equation 18) 
y’{t) is the radius of the MPD at time t. 

For any given moment t, 0 <t<T, u'{t) is the optimal constant speed before t to 
maximize the distance from the datum. The helicopter does not guess any particular value 
for t and makes the pessimistic assumption that for every moment in the interval [0,7] the 
submarine is able to reach y’(t ) . 



U’{t ) is an infeasible “Virtual Speed” function to cover every possible use of the 
submarine’s battery. 

Finally, the probability of detection when the helicopter searches within the MPD 
can be expressed as 



The following graphs show the evolution of some of the variables in Programs 1 
and 2 when considering the following realistic data: 

Constants: 




dt 



Then, 



U'{t) = -EroT n t 
3 



l/3,-l/3 



{Equation 19) 



a(x)dx 



PD ’ (/) -\-e 



y (x ) 2 



=;r <t<T 



{Equation 20) 



Helicopter’s delay t = 10 minutes. 
Sweep width W = 3 NM. 
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Submarine energy Etot= 5293 miles 11 hours' 1 . 

Helicopter’s speed V = 60 knots. 

6 = 0.05 hours. 

T= 2.83 hours. 

Optimized Variables: 

Initial distance submarine-datum y(r) = 4.2385 NM. 

7= 65.7236 NM. 

v = 5.64847 knots. 

Figure 2 shows the distance y(0 from the flaming datum to the FOC that 
minimizes the probability of detection for the submarine (Program 1), and the MPD 
considered in the helicopter’s search (Program 2). Euler’s method is used to solve 
Equation 14 numerically. The distance y(r) is utilized to compute the submarine speed 
u(t ) . Then, the computed speed is used to obtain the distance y(t + At) considering that: 
y(t + At) = y(t) + u{t)At, where At =0.02 hours. ( Equation 21) 



Distances FOC and MPD 




Figure 2. Submarine distance as a function of time. 
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Figure 3 shows the submarine’s velocity functions U\t), u(t) and u’(t). u(t) is 



the velocity function to reach the FOC shown in Figure 2 (Equation 14). U\t) is the 
submarine’s virtual speed that determines the MPD assumed by the helicopter (Equation 
19). Finally, for the function u"(t) every value is independently feasible, and each one 
implies the total consumption of the battery at time t (Equation 1 7). 



u(t), lT(t)andif(t) 




Figure 3. U\t), u (/) and u'(t ) when the helicopter’s delay is 
10 minutes and T= 2.82 hrs. 



Figure 4 shows the cumulative distribution function (CDF) for the time to 
detection as a function of time when the helicopter’s delay is 10 minutes and the search is 
executed: 

Inside the MPD circle (Equation 20). 

Inside the FOC (Equation 1 5). 
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Figure 4. Detection probability if helicopter’s delay is 10 minutes. 




Figure 5 shows the decrease in probability of detection as a function of the 
helicopter's delay x, if the search is executed inside either the FOC or the MPD. This 
graph is obtained from the analytical model for the parameters listed in page 15. 



Analytical PD vs Delay 




Figure 5. Analytical Probability of Detection versus Helo’s delay. 
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To obtain a more accurate result a discrete event simulation using the Java-based 
Simkit Software package is developed next, which considers the same characteristics for 
the submarine and helicopter used in the analytical model. The next chapter shows the 
methodology used to solve some of the problems related to the movement and strategies 
for the searcher and the target. 
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III. THE SIMULATION 



A. CONSIDERATIONS 

1. Submarine Strategy 

a) Submarine Initial Position 

The initial distance from the datum y{r) is an optimized variable for each 
value of x. To determine the initial position of the submarine for each replication, the 
distance y(j) is multiplied by a different right triangular random variate [0, 1] called to. 
The pair (y(r)cr, 0) is then transformed to rectangular coordinates and thus creates an 
initial position uniformly distributed inside the circle of radius y(r) . 



b) Velocity function 

The submarine’s general velocity function is generated through 

Equation 14 



11/ 3 



u(t) = V 



Y 2 



y(t) + 



2 , T(0 3 



-1 



; r <t<T 



SV 



where Y and v are optimized constants for the particular helicopter’s delay rand y(t) is 

the submarine’s distance from the datum at time t, obtained from Equations 14 and 21. A 
vector holds the velocity whose value is updated every 9 seconds in order to create a 
function which is relatively continuous. For each replication, the speed function is 
multiplied by the same to used to determine the submarine’s initial position. In doing so, 
a different speed behavior is generated for each simulation run and it is consistent with 
assumption a) for the submarine in Chapter II. 
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c) Course 

The submarine flees from the datum with a random course 9, restricted 
turn capability and executing 9-second legs. When a leg has finished, the submarine 
decides if it will change course. This decision is made by sampling a uniform [0,1] 
random variate and comparing it with a reference value. This reference value is very 
close to 1 and the maximum change in course is restricted to 3 degrees port or starboard. 
In doing so, the submarine maintains a fleeing attitude. 

2. Helicopter Strategy 

a) Position of the First Dip 

The helicopter’s first dip is uniformly distributed in the circular area 

whose radius is y'(r + time to dip), where y'{t) is the MPD function obtained from 
Equation 18. 

b) Position of the Next Dip 

The helicopter dips uniformly in the area whose radius at time t' is y'(t'). 
Washburn and Thomas [Ref. 6,p. 420] state that, considering X as the present helicopter’s 
position in the unit circle and t as the present time, the earliest time t ' to reach a point X 
in the same circle guaranteeing uniform search can be computed as: 
t'= t(d -abc + e)/(\ -b 2 ) = T(X,X',t) ( Equation 22) 

Where 

U\X\ , U\X'\ _ . . , „ . 

a = — — — b = — — — c = Cos {angle between X and X ) 

d = \ + — and e 2 - a 2 +b 2 d 2 -a 1 b 1 {\-c 2 )-2abcd 

t 
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U and V are the submarine and helicopter’s velocity and D is the time required for the 
helicopter to make a dip with the sonar. 

The result stated in Equation 22 makes it possible to compute the position 
of the helicopter’s next dip and to be certain that it will be uniformly distributed in the 
submarine's possible positions disk. Even though this procedure is applicable when 
dealing with a submarine with a constant maximum speed, a similar approach can be 
done when considering that for any time t , the submarine’s virtual speed value U' (r) 
may be obtained from Equation 1 9. 

Let 

X = [x,y] X'=[x',y'] 

be the helicopter’s position coordinates in the unit circle at time t and ( respectively, so 
x 2 + y 2 < 1 and x' 2 +y' 2 < 1 . 



Then, the time for the next dip t' should satisfy 



f , r (>•)-*>•• tor +wy\n -&(.<)? 

V 



( Equation 23) 



Where ^][x' y’ (t') - xy*(r)] 2 + (^') - >y ; '(0] 2 is the distance between the helicopter’s 

present dip position at time t and the future dip position at time Even though Equation 
18 allows obtaining the value of replacing it in the equation above creates a sixth 

degree equation that is not useful to solve for t'. However, we can obtain the approximate 
distance between [x, y] and [x', y'] if we assume that an approximation for y*(t') is: 



y'(0 = y(t) + irm'-t) 



where U’(t) is the virtual speed introduced in Equation 19. 
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Replacing y'(t') in Equation 23 we have: 



, + Dh + +[/(v-(0 + C/'(0«'-0)-x>- : (oF 

V 



(Equation 24) 

If both sides of the equation are multiplied by V, then 

vt+vD + Viy (/(o + £/•(/)(/’-/))- + ty (/(o + i/ w-o)- w (0] 2 = vt 

(Equation 25) 
and 

^ + d - n = -Vt4/(o + £/*(/)(/•-/))- */(oi 2 + [/ (>■ *(o + irm-t))- 

(Equation 26) 

Squaring and reducing, then 

V\t + D - t') 2 = [*’(/(') + i/ W-O)- */(0] 2 + b'(y*(0 + U'(t)(t'-t))- yy\t )] 2 

(Equation 27) 

Now let 



t + D = a, y(t)(x'-x)-x'U'(t)t = b, and y(t)(y'-y)-y'U\t)t = c. 
Then 



V 2 (a-t') 2 = [b + x'U'(t)t'] 2 +[c + y'U'(t)t'f (. Equation 28) 

and reducing 



i t a 



V 2 a 2 - 2aV 2 t'+V 2 t' 2 = b 2 + 2bx'U‘ (t)t'+(x'U‘ (t)) 2 1' 2 +c 2 + 2cy'U' (t)t'+(y'U' (t)) 2 t 



(Equation 29) 
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Now factorizing 




(Equation SO) 



If 

d = V 2 -(x'Lf(t)) 2 - (y'U* (t)) 2 , e = 2aV 2 +2bx'lT(t) + 2cy'ir(t) 



and / = V 2 a 2 -b 2 -c 2 



then, dt' 2 -et'+ f = 0 



(Equation 31) 



which is a normal second degree equation. Finally 



e±‘Je 2 ~-~4df 

2d 



(Equation 32) 



In practice we choose the smallest positive root of Equation 32, if greater than the 
current time t. The chosen time will be the solution for the earliest time that a point in the 
unit circle can be reached to guarantee a uniform dip location. 

Once t' is obtained , y*(t’)can be calculated using Equation 18. Therefore, the 
coordinates of the next helicopter’s dip are: 



Figure 6 shows two typical sequences of dip times when the helicopter’s delay is 
10 minutes. 



X'={x'y(t'), y'y(n ) 
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HELICOPTER DIP TIMES 




Figure 6. Helicopter’s typical dip times when delay is 10 minutes. 



3. Simulation Inputs 

a) Parameters 

In Equation 14, the parameters Y, v,and y{z) are optimized for different 
values of r using a spreadsheet and a non-linear programming solver to solve Program 1, 
as explained in Chapter II C. The following table shows the parameters introduced in the 
simulation for each value of r. All the columns except for the first were calculated based 
on the following constants: 

Helicopter dip radius R = 1 .5 NM. 

Submarine energy Etot= 5293 miles 1 , hours' 2 . 

Helicopter speed V= 60 knots. 

Time to make a dip 8= 0.05 hours. 

At for the spreadsheet = 0.02 hours. 
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Helicopter endurance = 2.65 hours. 



PARAMETERS INTRODUCED IN SIMULATION 


Delay 


Initial Distance 


V 


Y 










10 


4.238500882 


5.64847041 


65.72360358 


15 


5.677708987 


6.00920271 


73.40442537 


20 


6.981682236 


6.2818038 


79.96614511 


25 


8.1932953 


6.50380744 


85.87639327 


30 


9.340396291 


6.72857452 


90.75760417 


35 


10.43320542 


6.87398188 


96.13915272 


40 


11.48289288 


7.01625994 


101.0523998 


45 


12.49630944 


7.14253496 


105.7943132 


50 


13.47854091 


7.2559127 


110.3977112 


55 


14.43347615 


7.3589631 


114.88283 


60 


15.36430179 


7.45277939 


119.2756334 


65 


16.2685256 


7.51082216 


124.0679351 


70 


17.16295785 


7.61865811 


127.8289081 


75 


18.03435462 


7.69262138 


132.0131104 




able 1 . Simulation Parameters. 



4. Simulation Validation 

a) Average PD (PDa) 

The average probability of detection (PD) for a particular value of the 
helicopter’s delay r is obtained if it is considered that for a helicopter dip at time t, the 
probability of detection might be computed as: 

R~ 

PD(i' h dip) = 7 (Equation 33) 

/(O' 

Let N be the number of dips made by the helicopter, including the detecting dip (if any) 
but none beyond it. Then the probability of detection PD for any simulation replication 
is: 
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(Equation 34) 



( N \ 

PD = E\ Y J PD(i' h dip) 

V < = ! J 

The expected value is estimated by replication. Letting N k be the number of dips in the 



Ar* replication and t kt be the time of the z th dip in the k * replication, the estimate for Q 
replications is: 



Average PD a (t,Q) = 




(Equation 35) 



B. BASIC RESULTS 

10000 replications were run for 14 different values of r. The simulation 
simultaneously computed Monte Carlo estimates for the value of PD A . Figure 7 (next 
page) shows the evolution of the probability of detection and its 95 % confidence 
intervals versus ras estimated by the simulation compared with PD a, when the 
helicopter’s search is executed inside the MPD. It can be seen that the values of the 
estimates PDa remain inside the confidence intervals for all values of r, confirming the 
accuracy of the simulation results. In the next chapter the simulation is used to predict the 
probability of detection when the behavior of the submarine and the helicopter are varied, 
for helicopter’s delay of 10, 20 and 30 minutes. 

Tables with detailed simulation results used in this chapter are included in 
Appendix A. 
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PROBABILITY 



SIMULATION AND PDA VS DELAY ( MPD) 



e— SIMULATION 95% Cl PDA 




Figure 7. Probability of Detection obtained from simulation and estimate 
PDa versus helicopter's delay, when search is inside MPD. 
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IV. RESULTS AND ANALYSIS 



A. MODEL FITTING 



Figure 8 shows the probability of detection as a function of the helicopter’s delay 
r according to the results obtained from the simulation, when the search is executed 
inside the MPD. A fitted function is also shown and the corresponding regression 
analysis is detailed in Appendix B. That curve is valid for the parameters listed in chapter 
III. 
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B. TACTICAL IMPLICATIONS 

Different variations to the problems were explored. The following results were 
obtained from 10000 simulation replications. A large-sample 100(l-a)% confidence 
interval for a population proportion p, if n is the sample size is: 



When we don’t know the value of p the most conservative approach is to use p = .5 , to 
maximize the value of the expression under the root [Ref. 9, p.289]. Therefore a 95% 
confidence interval is: 



1. Submarine on the Edge of the FOC 

The submarine might decide to minimize the analytical probability of being 
detected by remaining on the edge of the FOC. Figure 9 shows the probability of 
detection for helicopter delays of 10, 20 and 30 minutes when the submarine escapes by 
sailing on the edge of the FOC compared with when its position is uniformly distributed 
in the same area. 

A small but statistically significant decrease in the probability of detection is 
observed for helicopter’s delay of 10 minutes (p value = 0.000975). 
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Slimarine cnFOC Edge Vs Uhiform 

□ Cn the Edge iLiTfcrm 




Helo Delay (rrirutes) 



Figure 9. Submarine escaping on the edge of the FOC compared 
with when its position is uniformly distributed in the same area. 



This effect seems to be attractive for the submarine. However, according to our model, 
this tactic has the disadvantage that at the end of the helicopter’s search period the 
submarine’s battery is completely exhausted and its options are reduced to being 
immobile, at the bottom or snorkeling. No statistically significant differences are 
observed for delays of 20 and 30 minutes (p values = 0.51 and 0.15 respectively). The 
tactical significance is small in all cases, which should not be surprising because the 
helicopter makes no attempt to exploit the predictability of the range to the submarine. 
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2 . 



Submarine Remains in the Vicinity of the Datum 



The submarine might attempt to surprise the enemy force by remaining in the 
vicinity of the datum. This situation may be also presented if the submarine executes an 
involuntary attack and is forced to remain close to the datum because of its low battery 
level. Figure 10 shows the effectiveness of this tactic for 3 different values of the 
helicopter's delay compared with when the submarine behavior determines a uniformly 
distributed position in the FOC. If the submarine remains in the datum vicinity, no 
statistically significant differences are observed for delays of 10, 20 or 30 minutes (p 
values = 0.31, 0.43 and 0.95 respectively). 



Submarine Remains in Datum Vicinity 

□ Vicinity ■ Uhiform 

02 
o.B 
0.15 
>* 0.14 

E 0.12 
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0.1391 



0.1342 




0.0713 0-0742 



El 



0.0503 0.0505 



10 20 30 

Helo Delay (Minutes) 



Figure 10. Probability of detection when the submarine remains 
in the datum vicinity, compared with when its position is uniform. 



These results indicate that the submarine does not get any tactical benefit 
remaining in the vicinity of the datum after the attack. Again, this should not be 
surprising given the helicopter’s tactics. 
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3. Helicopter Fails First Detection Success 

When the helicopter is dipping, a poorly trained crew may fail to detect the target 
even when the target is within range. When dealing with a diesel-electric submarine this 
is especially possible. They almost make no noise and most of them are small. Although 
it is logical to suppose that this lack of skills will negatively affect the probability of 
detection, our intention is to measure its real impact in the mission accomplishment. In 
Figure 1 1, a dramatic loss of probability of detection in the event that the helicopter fails 
to detect the first time the submarine is detectable, for all helicopter’s delays. 



Helo Fails First Detection Success 

□ detect second ■ detect first 

0.16 
0.14 
5 * 0.12 
= 0.1 
re 0.08 
O 0.06 
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0.02 
0 



Figure 11. Probability of Detection when helicopter fails to detect 
the submarine when an opportunity is presented. 

In this case it is very unlikely that the helicopter will ever have a second chance of 
detection. 



0.1342 




10 20 30 
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4. Helicopter Searches Inside the FOC 

The helicopter's search may be executed inside the circle determined by the FOC 
instead of the MPD. In that case, the helicopter assumes correctly that the submarine will 
remain inside the FOC the whole time and the search is conducted according to this 
assumption. When this assumption is correct, a small increase in probability of detection (p 
value s 0) is observed for all three delays (Figure 12). 



Searching Inside the MPD or FOC 

□ IVPD ■ FOC 

0.25 t ■ ■— 




10 20 30 

Delay (minutes) 



Figure 12. Probability of detection when the helicopter's search is inside 
the MPD or the FOC. 



However, assuming that the submarine is inside the FOC will be incorrect whenever the 
submarine exceeds the limits of that circle. To do that the submarine could follow any of 
many speed behaviors determined by «*(/) given by Equation 17. For these cases 
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neither the analytical model nor the simulation allow us to know the effects in terms of 
the probability of detection of searching based on a wrong assumption. 
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